NPS-53Ru73091 A 



NAVAL POSTGRADOATE SCHOOL 

Monterey, California 




PRELIMINARY RESULTS CONCERNING THE IMPROVEMENTS 
REALIZABLE THROUGH THE USE OF VARIABLE THRUST 
TOGETHER WITH ENGINE GIMBALING FOR A PARTICULAR 
INTERCEPTOR MISSILE 

by 

I. BERT RUSSAK 
SEPTEMBER 1973 

Approved for public release; distribution unlimited. 



FEDDOCS 

D 208.14/2:NPS-53RU73091A 




NAVAL POSTGRADUATE SCHOOL 
Monterey , California 



Rear Admiral Mason B. Freeman M. U. Clauser 

Superintendent Provost 



ABSTRACT: 

Some interceptor missiles as presently formulated possess a programmed 
thrust magnitude history with a gimbaled engine to provide steering. We 
examine one such missile to determine whether performance can be improved 
if we allow a variable thrust magnitude together with engine gimbaling to 
provide control. 

Two trajectory optimization programs were written to provide an initial 
answer to this problem. Preliminary results indicate reductions in the time 
to intercept by as much as thirty per-cent over that obtained by the pre- 
sently used guidance scheme. With tuning of the programs it seems rea- 
sonable to expect even greater improvements and further investigation seems 
warranted. 
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Introduction 



Some interceptor type missiles as presently formulated possess 
programmed thrust magnitude history with a gimbaled engine to provide 
steering. The present guidance scheme used on these missiles determines 
the steering control and hence the direction of the thrust vector. We 
examine one such missile and answer the question as to whether perfor- 
mance can be improved if we allow a variable thrust magnitude together 
with thrust direction to be controlled by some guidance scheme. 

In order to take the first step in answering this question, two 
trajectory optimization programs were written. These were designed to 
determine optimal histories of thrust magnitude and direction in order 
to obtain minimum time to interception for our missile under given 
scenarios. While the programs are not in a finely tuned state, neverthe- 
less, preliminary results indicate reductions in the time to intercept 
by as much as thirty per cent from that obtained by the present scheme. 

With tuning of the programs it seems reasonable to expect even greater 
improvements and further investigation seems warranted. 

Model 

The missile model used was two dimensional since all test trajectories 
were flown in a horizontal plane. 

Letting the indicated terms have the meaning specified in the 
nomenclature, then the picture of the model is: 
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Figure 1 
Missile Model 



4 



(1 N 

The differential equations for this model are : 
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where 


i) 


y x , . . . ,y 5 are called state variables since they define the state 



of the missile and TT,AT are called control variables since they control 
the state through the equations 1); ii) FA,FN are functions of the velocity 
vector and the control angle AT, 

The constraints for this problem are 
2a) 0 < TT < 14400.0 

TF 

2b) / TTdt < 38,500 

0 

in which 2a) is a thrust level constraint which says that our thrust must be 
non-negative and is bounded above by 14400 lbs. and 2b) is a condition on the 
amount of fuel used. 

Our task is, given the initial conditions 
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for the target, then determine a history of TT, AT in time which 
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Detailed equations are presented in the Appendix 
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yields a minimum for the time of intercept TF. Using the penalty method to 
include the constraint of target impact in the cost function, our cost 
function is then 

4) c = tf + un[ ( y 1 - y 1T > 2+ (y 3 - y 3T ) 2 ] 

Method of Solution 

A. General Techniques Available 

There are many ways to attack a problem of the type specified 
above. For example; 

a) the classical calculus of variations technique 

b) gradient technique 

c) conjugate gradient technique 

Of these a) is an indirect method, which seeks a trajectory which satisfies 
certain necessary conditions rather than seeking to reduce the cost function 
directly. This method depends upon the choice of the initial values of 
a set of multipliers called adjoint variables which satisfy a certain system 
of differential equations. This choice is often a highly sensitive one and 
instability in attempting to converge to a solution trajectory can result. 

Methods b) and c) are direct methods in that they directly seek to 
minimize the cost function by seeking new trajectories with lower values 
of cost function. All of these methods are based on generating a sequence 
of trajectories which converges to the minimizing one. The gradient 
technique works by linearizing the cost function at each trajectory of 
the sequence developed and iterates to the next trajectory of the sequence 
by changing the controls in the direction opposite to the gradient. The 
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conjugate gradient technique is a step more sophisticated than the gradient 
technique in that it generates new trajectories in its sequence by effectively 
expanding the cost function in a Taylor Series up through the second order, 
thus obtaining a more accurate representation of this function. 

All of these methods together with a number of others were considered 
for the problem at hand and because of greater sureness of convergence 
the conjugate gradient method was selected. 

B. Brief Description of the Conjugate Gradient Technique 

This method is most easily described when discussing the problem 

of minimizing a cost function which is a quadratic function of the N variables 

x. , » • • »x . Thus assume that we are given the problem selecting values of 
1 n 

x^,...,x n in order to obtain a minimum of the quadratic function 
5) c(X) = d + BX + 1/2 X T GX 

where: i) X denotes the vector (x-,...,x ); ii) d denotes a constant and B 

In 

denotes a constant vector; iii) G denotes the matrix of second partial 
derivatives of c. Given a starting point the conjugate gradient method 
computes a sequence of vectors Hq, H^,..., along which the function c is 
minimized. Thus starting at Xq the method computes a direction which 
depends on the cost function c and the point X^ and determines a value X^ 
which is a minimum of c in that direction. Next, a direction is 

computed at X^ and c is minimized along that direction to produce the point 
The sequence continues in this manner and it can be shown that in the 
absence of round-off, the method will converge to the minimum point in at 
most N iterations (where N is the dimension of the vector X). 
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In general, as in our case, the cost function is not quadratic. The 
procedure then is to approximate the cost function by the first three terms 
if its Taylor Series at each iteration point so that it has the form of a 
quadratic and to develop the directions from those approximations as 
outlined above for the quadratic case. Details of the conjugate gradient 
method as originally developed by Hestenes for linear systems, are in [1] 
and its application to general functions is explained in [2]. Furthermore, 
the technique of conjugate gradients works on more general functions than 
functions of a finite number of variables and one may apply it with some 
modification to functions of an infinite number of variables (see [3]). 

Thus for a cost function which depends upon an infinite number of variables 
as our cost function which depends upon the value of TT and AT at each 
time point, one may use this technique to seek out those values which 
minimize it. 



C. Application of the Conjugate Gradient Technique to Our Problem 

In order to apply the conjugate gradient technique to our problem, 
two computer programs were written. 

The first of these programs was written using the conjugate gradient 
technique for an infinite number of variables as referred to above. This 
program is listed in the Appendix B and was never fully checked out due to 
lack of time. 

The second program was written using the conjugate gradient method 
for functions of a finite number of variables as outlined above. Now as 
previously stated, the cost function for our problem depends upon infinite 
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dimensional controls, namely the magnitude TT and direction AT of the thrust 
vector at each time point • However in any computing machine procedure f or 
integrating the differential equations for our problem, only values of 
the controls TT and AT at a finite number of time points are used. For 
example, in the simplest type of integration scheme, if the time interval 
is denoted by DT and t Q , t^, • • *ty . . are the time points of the 

integration scheme then 

y(t x ) = y(t 0 ) + y(t 0 )-DT 
6) y(t 2 ) = y(t 1 ) + yCt^'DT 

y(t +1 ) - y(t ) + y(t )-DT 

J J J 

y(TF) = y(TF-DT) + y (TF-DT) *DT 

where y,y denote the state variable to be integrated and its derivative and 
TF denotes the final time. In this scheme only the values of TT and AT at 
the time points t^ affect the trajectory. Thus our cost function which 
depends upon y at the final time in turn also depends on the values of TT 
and AT only at these time points. 

Thus, the computer really reduces the infinite dimensional problem 
to a finite dimensional one. Furthermore if we take this into account 
in formulating our model then our numerical optimization scheme which must 
abide by such shortcomings of the computer, will be surer of success. 

This then is the technique used to adapt the finite dimensional 
conjugate gradient method to our problem. The integration scheme selected 
is the one used on already existing trajectory computer programs for the 
missile under consideration and is as follows: 
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y 5^ t j+l^ * f 5^j>' DT 



where we have denoted by f i = 1,...,5 the right hand sides of 1) . This 



by using the first two derivatives of position, while integrating the 



vatives of these quantities. 

Besides computation of the cost function at each iteration point, the 

conjugate gradient method requires us also to compute the derivative of 

the cost function with respect to the control variables TT(t^), AT(t^). 

By the chain rule for differentiation, this requires that we first 

differentiate the cost with respect to the state variables at TF and then 

differentiate the state variables at TF with respect to the controls at 

the times t.. The former derivatives are easily formed, however the latter 
J 

derivatives are formed sequentially as follows: According to the inte- 

gration scheme 7) forming the derivative of y^(i = 1,...,5) at tg with 
respect to AT(tg) and TT(tg) yields 



integration scheme essentially integrates the position components y and y^ 



velocity components y^, y^ and the mass y^ by using only the first deri- 
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Forming the derivative of at 
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with respect to AT(t^) and TT(t^) yields 
STTCtj^) ° 1 



and next, forming derivatives with respect to AT(t Q ) yields 
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i = 1,3 



where the last equalities in 9b) result from 9a) and where f.(t. ) means the 

1 ic 

function f^ evaluated with arguments y(t^), AT(t^), TT(t^). Similar equations 
hold for the derivatives with respect to TT(t^) and TT(t^). Continuing in 
this fashion, then at time t^ we form the derivatives of with respect 

to TT and AT at all time points up through t^. Forming the derivatives 
with respect to AT at all such times, first we set, (as in 9) the derivative 
with respect to AT(t^) 
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and finally, for the derivative with respect to AT at all other preceeding 
time points t gl s = 0,1,.., k - 2 
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with similar equations holding for the derivative with respect to TT at 
all time points* It is recognized that all derivatives of the state y 
required on the right hand side of 10) have already been formed at previous 
steps in the process* 

This procedure continues until we reach TF and thus obtain the 
required derivatives of final state* 
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Since the cost function also depends upon TF, we are also required to 
form the derivative of the cost with respect to TF, however this presents 
no difficulty. 

Results 

In stating the results of using the above described computer program, 
it is to be noted that there were severe time limitations on this initial 
phase of the project so that only a minimal amount of time was left after 
formulation, development and checkout of the basic computer program. 
Consequently, the results presented herein are preliminary in the sense 
that no "tuning" (such as problem scaling) of the computer program to 
this problem was done. Such tuning will produce better and often very 
significantly better results than the basic program. Nevertheless, the 
results that were obtained indicate significant savings in time over those 
obtained from the presently used guidance scheme. 

The basic missile target scenario that was used had the target at 
20,000 feet initial range. Both missile and target had initial velocity 
of 800 feet per second. The missile heading and target aspect were varied 
as depicted by dashed lines in the figure below. 
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Figure 2 

Basic Missile-Target Scenario 



The results of the conjugate gradient runs together with a comparison 
to the results of the presently used guidance scheme' are presented in the 
table which follows. In addition, plots of some of these comparison 
trajectories are also presented. In each plot is indicated the time of 
intercept with the target. Finally, the values of the control variables 
TT and AT at each time point t^ are listed for each plot. The number of 
such time points or equivalently the number of intervals in the integration 
process is arbitrary and was generally selected to give roughly an interval 
of .25 sec for the initial trajectory and time of flight which were used 
to start the program for each case. 
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Comparison of Times to Intercept Obtained By Conjugate Gradient and Presently Used Scheme 
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History of Thrust Magnitude (lbs.) and 
Direction (Radians) at Each Time Point 



Missile Heading -90° 
Target Aspect 90° 



THFMJST USED A»M6LK USED 

-A . 605^09 5*>1 4 595992-06 A .733 1 2 1 593 v . 07<'7 

1 4,2 9 9 . 9 9 9 9 7 9 3 3 v 9 .30598 4 3 9 0 3 5 8 4 4 3 3 

14399.99997999224 .3803681 1 171 61034 



1 /■ 3 9 9 . 9 9 9 9 8 0088 y 6 
1 4399 • 99 9980 3 P 638 
1 4399 . 9999? 06 08c 6 
14399.999980874 
14399.999981 17321 
14399.99998147072 
14399. 99998 178906 
14399.0 9 9 0 88. 1 8 8 6 ; v 
14309 .9 8? .•."284 
1 4399 .99998298058 
14399 .99998352771 
1 4399 .9999841 7795 



•376756521 4444642 
•3748835541 679647 
.374 2 06295 5 >'• ? 6 3 2 .s 
.3739137777928477 

.371 6230858999721 
.3662140979628795 
.361 101 9836883072 
•356230888 1 0 1 iOOv 
.3 51 6 1 4/122 7(131 92- 
.3 4722 330 46 08 1 528 
. 34298 7 42 51066 72. 6 
.33879451 48 1 94222 



14399 .9999849742 
-1 . 48 4 1 3 5400 174 1 98E-0 5 
-1 .367361742331 171 E-OS 
-1 .29380 1 1 •• 5470 2 59 5E- 0 5 
-1 .201 62303 1 57 587 4 £-0 5 
-1 .112 5974 41 37462 yr>0 5 
-1 .9261 i t i - i •• • /■-; •:■ i -95 
-9.42 1 1 664930391 558-06 
-8 • 6 0 6 2 4 7 0 8 3 2 3 3 2 3 5 E - 0 6 



3344809309495753 
.3395337541740675 
*3383031 o' l 5055445 

• 33747 43 7 2 9 38573 9 
*336939 1 8 8487 191 

• 33 6 ' 4 5 v 9 >401 62 '9 
.' .o' /.9701 01 92 663 

'♦33649044041 70523 
•3366156214075165 



-7.8161 166 53 9 9 1 8 0 6 E-0 6 
-7 *0506530 59250 1 09E-0 6 
-6 *3097 543391 440738-06 
-5 . 59333 7347 902 79 7F -06 
-4.90133603^8383 1 5E-06 
-4 *23370489874 50 38 E-06 



.336 8, 4 9 7 5 1- 5 6 1 40 6 4 
,..33718591 12 59483 
.337021 6675022808 
.3381 581899218649 
.3307996 66 1726958 
.3395529221752953 



-3 . 590 41071 386957 4E-0 6 



.3404271763094834 



-2 .97 1 4405670 82 99 6E-0 6 
-2 .376798128100541 E-0 6 
-1 * 806 5050 1 1 77/- 5 49 E-0 6 
-1 .26060166 140531 4E-06 
-7 .39 1/1863 539001 89E-07 
\ -2 .422264241 52342 1 1>07 
V v 0 •35000000000 ci 3 6 - 8 



• 341 4338 744 12035 
.342566 560 4 4 7 3136 
.34390 O' 7 3 7 8455 /•• 1 8 

• 3453936 690 1 On 5 7 
.34708404 / 0 9 0 0013 

• 3489v 142 1 94-401.1 
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History of Thrust Magnitude (lbs.) and 
Direction (Radians) at Each Time Point 



Missile Heading -90° 
Target Aspect 180° 



1:1: AST USED 40GLE 
-1 • 1 646*6*1 5P55909E-C 
1 A3 9 9 .99998027 ni9 
1 4399 .9999^0475*4 
1 4399 .99993069? i /. 

1 4399 .99993079 94 
1 4399 »9 9 9 9 9 0 s' 0 4 67 
1 439 9 *9 9 j 9 0 ' : 9 5 6 1 6 
14399 .999931 1243 5 
1 4399 .99 993 1 32 3 58 
1 4399 .99993155864 
14399.99998185297 
1 4399.99993222639 
14399.99998270991 
14399.99998334043 
14399.99998416803 
-1 .566493161 176776E-05 
-1 . 46087493 5475238 £-0 5 
-1 .357041 566137291 E-05 
-1 .256871 182486963E-05 
-1 . 1 598228294433 53 E-0 5 
-1 .0659421870431 55E-05 
**9.7 52 6 57 448003 036E -06 
-8 .87323856391 425? r -06 
-8. 03 64300220 162 77 E-06 
-7 .227474929841 107E-06 
-6. 45 1 60777302 1927E-06 
-5.7090646P9795342E-06 
-5. 000091 267843396E-06 
-4 .32 494 9 57 73 907 57 F. -06 
-3.683922571 863313E-06 
-3 .07 7 04360986 82 48 E-06 
-2 .504968690771 787E-06 
-1 .968061221 1 1 89 1 7 E-0 6 
-1 .4701 7424471 0638E-06 
-1 *007 1 27 876 053239 E-0 6 
-5. 7932633960681 6E-07 
-1 • 871986448 563937E-07 
0 0 



USED 

5 4.723242003744533 

6.47 9 2 6 2 2 9 4 9 ? 7 5 9 7 E - 0 2 

6 .33 528543 69 506 43 E-02 
6 .437 6 9 0 9 6 1 8 6 7 1 9 E - 0 2 
6.7971362091 13957E-02 

7 .4] 51697 520499 5 5E- 0 2 
7 .622 1 22 7 00236574 E - 0 2 
7.7008763662 12161 E-02 

7 *8285094882 961 9 3 E-02 

8 .01 4004523567484E-02 
.27357 58 5592 3 544E-02 
.6328254731 8824SE-02 
.13 1522395380891 E-02 
.”832592 52 5874304E-02 
1083958274495448 

7 .969 501 83 1 33 4 1 9E-02 
7 . 59503 74 1 96 7636E-02 
7.3014551 4 5550 6 73 E-02 
7 . 06 79569 74.3630 7E-02 
6. Q 781 89976333292 E-02 
6.7191 4507 7868842 E-02 
6 .58020 1445930 51 6E-02 
6.4523159517 9 40 5 7 £-02 
6 .32743446844/ 71 8E-02 
6.1 979758741 7582 7E-02 
6 *0564 1 2099 5 5062 4E-02 
5 .8948921 71 735704E-02 
5.70489057 1965851 E-02 
5.47686 5053 7802 68 E-0 2 
5.19991 707321 9 558E-02 
4. 86 97 94646 4643 6 5E-02 
4 . 4440980 55 129932 E-02 
3. 9256 4 571099 52 7 E-0 2 
3 .2726541 422 16548L-02 
2 .621487103^441 77 E-02 
1 .76851 454136273BE-02 
6» 5133138 59 402577L-03 
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